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Abstract. We calculate correlation functions of the (signed) density of zeros of 
Gaussian distributed vector fields. We are able to express correlation functions 
of arbitrary order through the curvature tensor of a certain abstract Riemann 
Gartan or Riemannian manifold. As an application, we discuss one and two- 
point functions. The zeros of a two-dimensional Gaussian vector field model 
the distribution of topological defects in the high temperature phase of two- 
dimensional systems with orientational degrees of freedom, such as superfluid 
films, thin superconductors and liquid crystals. 



PACS numbers: 02.40, 47.32, 42.30 

1. Introduction 

Topological defects play an important role in two-dimensional physics, and are 
observed in liquid crystal films 1 , two-dimensional crystals , superconducting 
films I^EIE], films of superfluid helium |7j, arrays of Josephson contacts 0111, nodal 
points of quantum and microwave billiards IllL [T^ and nodal points in optical 
waves ^1 ^1 . Common to these systems are orientational degrees of freedom with 
a global 0{2) symmetry. The above systems are either in an orientationally ordered 
or disordered state. Even in the ordered state, the orientational order is usually not 
perfect, but modified by the presence of massless spin waves. The corresponding 
correlation functions decay algebraically with distance (quasi long range order). In 
the disordered state, orientational order is destroyed through the spontaneous creation 
and subsequent unbinding of pairs of topological defects. The disordered phase is 
characterized by a finite density of these defects and exponentially decaying correlation 
functions. A suitable order parameter for two-dimensional orientational order is a two 
component vector field u{r) = (ui, ■ti2)(ri, r2). Topological defects are points, where 
the amplitude of the director field u vanishes. Circling around a defect, the phase 
of u adds up ±2tt (or rarely multiples of ±27r). We distinguish between positive and 
negative zeros (defects), where the sign of the defects (sign of the phase jump) is equal 
to the sign of the Jacobian det{diUj) right at the zero tt = 0. Figure ^ illustrates 
the two types of zeros. A crucial feature of the defects is the topological constraint on 
the number of positive zeros minus the number of negative zeros on a closed surface 

m 

number of positive zeros - number of negative zeros = 2(1 — 7) (1) 
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Figure 1. The planar vector field (ux 
(left) and a positive zero (right). 



Uy) = (x^ — 1,1/) with a negative zero 



where 7 is the number of handles of the surface. A sphere has therefore at least two 
positive defects, whereas a director field on a torus might be free of defects. Additional 
defects are created in pairs of opposite sign in order to obey the topological constraint 
In this paper we consider the distribution of defects in the disordered (high 
temperature) phase, where the director field u has a Gaussian distribution. At first 
we calculate general correlation functions of the (signed) defect density for arbitrary 
vector fields with a Gaussian distribution. As an application we study the mean defect 
density for certain distributions of u on curved manifolds. Then we give a transparent 
derivation of the two-point correlation function for a particular class of vector fields. 
Finally we check the neutrality of the defect distribution. 

Densities of zeros were discussed before in the context of nodal points of 
chaotic wavefunctions in microwave billiards [Ifil 1171 I12j . random surfaces ^1 
m EDI El El EHl El E3 EHl E3 and for the XY-model EHl Ei- Extensively 
studied are the optical dislocations (optical vortices) of "chaotic" optical wave fields 
IT3[TlinniiniE3lM,^i^^.37»i^ElEnillIl which are weU modelled by Gaussian 
distributed wave fields. The distribution of zeros of Gaussian random polynomials and 
Gaussian random analytic functions was discussed in |42[ 1431 1441 1451 146L WI\ . Random 
polynomials and random analytic functions have only positive zeros| and are therefore 
quite different from the systems considered here. 

2. Correlation functions 



We denote the positions of the zeros of a d-dimensional vector field Ui in d-dimensional 
space by Va and the charge (type of the extremum) by qa = signdet {diUj{ra))- The 
charge density reads p{r) = QaS'^ir ~ fa) and can be written as 

p{r) ^ det id,u,{r)) S^'iuir)). (2) 

To prove this representation, consider a particular field Ui with a zero at the origin 
and its expansion Ui{r) = Aijrj + . . . for small r, where Aij = djUi{r = 0) is the 
derivative of Ui at the origin and summation over double indices is implied. Then 
dei{diUj) = det A + . . . and 5'^{u) = S'^{Aijr.i) | det A\^^S'^{r). Combining both, we 
obtain 

demu,)6''iu) = ^^^S'^ir) = ±'^'(r). (3) 

J The real and imaginary part of an analytic function are harmonic functions, which have only saddle 
points as extremal points. 
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We consider now a Gaussian distributed vector field with correlations Xij 

(r,r') = 

{ui(r)uj{r')) and provide a scheme to calculate arbitrary correlation functions 
(p(ri) . . . p{rf)) of the charge density p, where ri, . . . , are / points and (. . .) denotes 
the average with respect to the Gaussian distribution of the field u;. 

We write a determinant as an integral over Grassmann variables 9i , rji 

detA = J{d7]d9fexp{9ir]jAij) (4) 

(for an introduction to the calculus of Grassmann variables see 03). To regularize 
the theory it is helpful to replace the sharp delta function by a Gaussian with a small 
but finite width uj and to use its Fourier representation 

SU^) = exp (-|!) . (2.)- / d'^pexp (-i.p^ - • .) . (5) 

We obtain 

lim(27r)"^ / d'^p{df]def exp {-up"^ /2 + 'i'{r,p,0,r])) (6) 



where the field 4* depends on 2d bosonic coordinates (r, p) and 2d fermionic variables 
{6, 77) and reads 

*(r,p, e, 7]) = ipiU,{r) + e,rijd,Uj{r). (7) 

Since ^ is linear in Ui and its derivatives, it has also a Gaussian distribution with 
correlations 

{^{xa)^{xb)) = {{efrjfd^u.irA) + ipfu^irA)) {O^ duui{r b) + ip'^uuirs))) (8) 
where x = {r,p,9,ri). The correlation function for finite width lo translates into 



(9) 



(p(ri) . . . p{rf)) = {2Tr)-f'^ J d^'^pidrjdey'^ 

X /expl-(^/2)^p2 ^^vl/(a;^) 

\ \ a a 



The Gaussian average can now be performed with the help of 



^exp ^{xa.)j ^ = exp ( i ^ (^'(a;„)*(.T^)) 



(10) 



yielding (summation over double i,j,k,l,m,n is implied, the partial derivatives act 
on the leftmost field only) 

{p{ri)...p{rf)) 

= I {drjdey'{2n)-f'' I df'pexp (^^J^'^i'^Pkri? {dr^MdkUiirp)) 

-Ij^PiPk M*fc<5a/3 + {Ui{ra)uk{rf!))) j . (11) 



a/3 



2 

a/3 
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After a straightforward (Gaussian) integration over the variables p we find 
(p(ri)...p(r/)) = (2^)-/'*/2(detg)-i/2 / (d,/d0)/'^ 




- {d,u,{r^)urn{r,,)) g™'''"'' {un{r,)dkUiirp))) j (12) 
fj.v / 
where we have introduced the formal metric tensor 

9ia,j/3iri, . . . ,r/) = LoSrjSap + {Ui{ra)u.j{ri3)) (13) 

and its inverse g'^°'-'^P — (gia.jp)^^ , where a pair of a roman and a greek index ia 
is understood as a single composite index. Now we interpret the tensor giajp as the 
metric tensor of a particular, / x c?-dimensional manifold. The manifold is parametrized 
by / X d internal coordinates a — {ri, . . . ,r f). It becomes singular in the limit uj ^ 
whenever at least two points Ta approach, since for two coinciding points two rows of 
the metric tensor become equal and det{gia,jf3) = therefore. These singular points 
show up as delta-like singularities in the correlation functions as shown in section 
In addition to the metric tensor we introduce an affine connection 

r -^..a = Sep {d.u,{rfs)um{r,)) 3"^''=^. (14) 
The affine connection allows us to define a covariant derivative of a covariant vector 
field Vjfj 

D,aV,p = 9„VS-/3 - T';i^^Vk^ (15) 
where dia = 9/9r''". The covariant derivative of a rank two tensor is§ 

Diatj/S^kj = diatjp^kj — ^ji3^^ia'^mfj.,k-i — ^ k-^^ia^j f3 ,mti ■ (16) 

A simple calculation shows the compatibility of the covariant derivative with the 
metric tensor Diagjf3,k-f — 0. The metric tensor together with the affine connection 
define a so-called Riemann-Cartan manifold, which is a generalization of a Riemannian 
manifold j^. Unlike the latter, it has a non-zero torsion tensor T^^jp — ^ia jp ~ 
■ The commutator of two covariant derivatives defines the curvature tensor 

(DiaDjp — DjfjDia) Vk-y = Rk-'/^iaJpVlx + pDmfj.Vkj ■ (17) 

A straightforward calculation yields 

Rkj,ix,ia,ji3 = SajS/3x {dtUk{r^)djUi{r x)) 

- SajS/sx (9iUfe(r^)Mm(r^)) g™'"'"'' {un{r^)djUi{rx)) 

- i^a) ^ UP) (18) 
where have used the metric tensor to transform a contravariant (upper) index into a 
covariant (lower) index giajpV^^ = Via- The curvature tensor allows us to simplify 
the expression for the correlation function 

(p(ri)...p(r/)) 

= (2^)-/'^/2(dct.g)-i/2y' (d,yd0)/'^exp (^-\9l9t^rv^ Rk^,ix,a.jp 
= (2^)-/'^/2(det.9)i/2y'(d,yd0)/'^exp (^-le^vtv^ R"^ ix'"" (19) 

§ The covariant derivative of higher-rank tensors is defined analogously, see 1151 . 
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It is quite remarkable that the correlation function is a purely geometrical object 
from the point of view of differential geometry — it depends on the metric tensor 
and the curvature tensor of a Riemann-Cartan manifold with metric tensor gia,jf3 = 
toSijSai3 + {ui{ra)uj{ri3)) and affine connection V'^^^^ = Sap {diUj{ri3)um{rp,)) ^'"^■'^t. 
It should be noted that correlations of the number density |17[ BUI (absolute value of 
the density) \p\ cannot be obtained by the above formalism. 

In case of a vector field Ui placed on a curved Riemannian manifold with 
metric Gij{a), one has to refer to a covariant version of the density of extrema 
(cr = (cr^, . . . , cr'*) are the internal coordinates to parametrize the manifold, is 
the covariant derivative related to the metric G^ ) 

p(cr) = det(G'^VjUfe)^cxp {^-^'^'^^i^J^ 

(det G)"^/2 det{d.,Uj)6^{u,) for cj ^ (20) 

where S/iUj = diUj at the zeros of Ui. The density differs from the Euclidean version 
only by the factor (detG)~^/^. 

2.1. Gradient fields 

A particular important class of vector field are the gradient fields Ui = di4> with a 
Gaussian distributed scalar field 0. One might think of as a Gaussian random 
surface |18l fl^' with local height z = (j){x, y). The zeros of 9^0 are the extremal points 
of (j). Extremal points with a positive signature are the maxima and minima of the 
function, whereas the saddle points of (j) have a negative sign. The scalar field <j) 
should not be confused with the (polar) angle of a vector field u = m (cos0, sin0). 
The angular field 9 is defined modulo 2tt and has singularities (at the zeros of it), 
whereas the field (j) is nonsingular and single valued. The metric tensor 113|l for di4> 
reads 

gia,jl3 = ^Sij^aP + {di(j){ra)dj(j){rp)) . (21) 

(The latter metric tensor is also introduced in [251 126L HT] in the context of excursion 
sets). The corresponding affine connection 

r -^,.0 = {drd,4>{rp)d,r.4>{r,)) g^^^^-* (22) 

is symmetric in the lower indices, i.e. the torsion tensor zero. The vanishing 

torsion characterizes a Riemannian manifold, where the affine connection is determined 
solely by the metric tensor g 

^ jh.ia ~ 2 ('^jl39mfi,ia + diagmfijp ~ dmfj.gia,j js) ■ (23) 

The Riemannian curvature tensor reads in this case 

Rk-y.lX.ia.jP = S/SxSa-y {{djdl(l){r p)didk(j)[r a)) 

- {djdi<l){r0)d^(j){r^)) g"'^^"'' {dn(p{r,)dA4>{ro.))) 

- (*a) ^ {j(3). (24) 

The defect density correlation function reads {p{ri) . . . p{rf)) ~ (27r)~-'^''/-^(det 5)^/^r, 
where the Grassmann integral 



r 



(dr,d0)/'^exp (^-i^^^fryf^^^E'^'r^V^^g) (25) 
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is the so-called total curvature. It is zero for odd dimensions (here: dimension — f xd). 
Its (covariant) integral over a (closed) surface is a topological invariant according to 
the theorem of Chern (1944) 0^ ^]. Indeed, the /-fold spatial integral over the 
manifold yields 

/ r 

n / d'^r„(det G(r„))i/2 {p{r,) . . . p{r^)) ^ (2(1 - 7))^ • (26) 

a=l 

We will now specialize our main resuh ijini. At first we calculate one-point functions 
and two-point functions. 



3. Density of defects the one point function 

3.1. The high temperature phase of the XY-model 

We calculate the defect density of a XY-model, placed on an arbitrary curved 
surface, in the disordered (high-temperature) region Pj . The Gaussian weight of 
a configuration Ui in the high-temperature region is 

P[u] cx exp ^-i j AA (V,Mj V*u^ + (r -I- (77 + l)i?/2) u^)^ . (27) 

where dA = -^/det Gijd^a is the invariant area element, Gy is the metric and Vi is 
the covariant derivative of the surface (a short introduction to the differential 
geometry of surfaces can also be found in (Sill)- The mass t determines the 
correlation length ^ = t~^/^ of the vector field and, therefore, the distance from 
the critical point. In addition, we have introduced a curvature dependent mass term 
as discussed in [211 ■ The weight H27|l is invariant under global 0(2)— transformations 
Ui —> cos(a)ui + sui{a)ei^ Uj , where eij is the covariant antisymmetric unit tensor with 
respect to the surface. The metric tensor reads gij = (uiUj) = dj {u,nu"^) /2, since 
Gij is the only symmetric rank two tensor invariant under 0(2). The affine connection 
is 

r^^, = {d,u,u„,) g"^" = (V,u,M,„) g"'^ + 7^, (28) 

where 7*^^ is the affine connection of the underlying surface. Due to the 0(2) invariance 
we can write H3 



(ViUjW„i) = ^Gj,n {ViUnU"^) + ^ej^e^P {ViUnUp) (29) 



and 



r,', = 7ji + e/f^. + Sp. log {{u')) /2 (30) 

where il^ — e"P (Viit„Up) / (u^). The gradient term does not contribute to the 
curvature tensor since it can be eliminated by a gauge transformation. The curvature 
tensor reads is the Riemann curvature tensor of the surface, R is the scalar 
curvature of the surface) 

Rkij — R^kij — (yi^j — V j^i)tk (31) 

and 

(^') 

Rijkl — „ — e™"Vmr2„) tM^ij- (32) 
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The Grassmann integral (|37|l is now trivial, and we obtain finally 

27rp = (i?/2 - e^^ V,%) . (33) 
The topological constraint for the total charge is easily verified since 

27r j dAp = j dA{R/2 - e'^V.,flj) = J dAR/2 = 47r(l - 7) (34) 

due to the theorem of Stokes [^0]. This expression for the density of defects was 
already obtained by a direct calculation and further analysed in [29) . It was shown 
that the signed density of zeros of the Gaussian XY-model with distribution (|?7|) is 
equal (up to terms of order t~^) to the charge density of the two-dimensional Coulomb 
gas in the Debye-Hiickel approximation. This is not surprising, since both approaches 
can be used to describe the high temperature region of the two-dimensional Coulomb 
gas In the next section we show, that the extremal points of a Gaussian field (j) 
behave in a similar manner. 



3. 2. Incompressible and irrotational vector fields 

We calculate now the defect density (p) for vector fields which are subject to an 
additional constraint, namely, the divergence or curl of the field must vanish. A 
field with vanishing divergence (curl) can be expressed via a 'potential' through 
Ui = dj<j) or Ui = dicj) respectively. We plug this representation of Ui into the weight 
H27|l and obtain the statistical weight for the scalar field 



P[0]aexp(^-i J dA{{A 



LB'/')^ 



r;i?/2)9,09» (35) 



where Alb = G'^^Widj is the covariant Laplace-Beltrami operator of the surface. The 
defects are the points where the vector field Ui is zero, i.e. where the potential (j) has 
an extremal point. The charge of the defect is the sign of the derivate of Ui: 

signdet(9iMj) = signdet(9i9j0) (36) 

for both cases. The mean defect density is according to theorem (|19|l 

p= {2n)-\dctGy^/^{detg)-^/^ J dTjideidr]2d92 

= MdetCy/^ j d'7ideidr/2d02 exp (^-i^fe^/ry.r;, R\^^ (37) 

where gij = {di(j)dj(f>). The Riemannian curvature tensor in two-dimensions has only 
one independent component-the scalar curvature R = i?*^' ki 

RkUj = {didk(f>djdi(f>) - (9,afe0a„(/)) g'"" (a„09j9/(/)) - (i ^ j) 

= ekie,jR/2 (38) 

where is the covariant, antisymmetric unit tensor with respect to the metric g^: 
ei2 = \/det g, £21 — —^/detg. The Grassmann integral H37|l is trivial, and we obtain 
finally 

27rp\/det G = {R/2)y/deig. (39) 
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The density of defects is according to equation H39|) the scalar curvature of an abstract 
surface (Riemannian manifold) with metric gij — {di4>dj(j)) . Therefore, with the help 
of the Gauss-Bonnet theorem and of equation (|39|l 



27r(number of positive defects — number of negative defects) = 2n J d^a\/ det 



9P 



d^ay/detgR/2 47r(l - 7) (40) 

where 7 is the genus (number of handles) of both the abstract surface gij\\ and the 
underlying surface Gij, since the abstract surface is closely related to the underlying 
surface and especially inherits its topology. The calculation of the curvature can be 
done in a large r-expansion (high-temperature expansion) for the special case rj = 0. 
From general considerations, the first two terms of the l/r-expansion of gij read 

ffy = = AG,, + B-G,j + 0(t-2) (41) 

r 

with constants A, B to be determined. We multiply equation 141() with G"-' and obtain 
(i? is the scalar curvature of the surface; the Laplace-Beltrami operator Alb refers to 
the surface as well) 

G^'g,, =2A + 2BR/T = (a>9,:<^) = ^ A^b (<^') - m.B<l>) 

= ;^Alb (((/>(t - Ai.b)0) + (0Alb0)) - (0Alb0) 

Zt 

= ^log(^)--^ + 0{r-') (42) 
An \a^T / 127rr 

where a is a short-distance cutoff. In 29 the following relations are reported 
Alb (0(r - Alb)(/') = - i?/(47r) 

]_ _R_ 
47r Va^T/ 247rT 
The density p can now be obtained using the fact that the metric gij is conformally 
equivalent to the underlying metric Gij (up to order 1/r). With the further help of 
equation H39(l we find as our final result^ 

^'^^-^/2 4^-l°<8l^°<i)-2l7) 

- 1 

2" + 6rlog(l/(aV)) 

R ^ . R / . . s 

-2+6^^-2 
where Z = log(l/(aV))/(27r). 

Expression (|44|l demonstrates, that the 'gas' of zeros of an incompressible or 
irrotational Gaussian vector field or likewise the set of extremal points of the potential 
(j) (random surface) behaves approximately as the two-dimensional Coulomb gas in 
the high temperature phase. The latter system is well described by the Debye Hiickel 
Hamiltonian for the continuous charge density p 



(0Alb0) =TZ^og[^]+^_+0{r-^). (43) 



^^^JdAjdA' i2np - R/2), G(a, a') {2np - R/2),. + ^ / - 

II The number of handles is zero for a sphere and one for a toroidal surface. 
^ The terms of order r" and do not depend on the coefficient 77. 



(45) 
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where x is the fugacity of the charges. By setting 5H/5p = 0, we obtain for the mean 
charge density 

^ ,_\^^ ^ - K + 4^^^ + (46) 

which is equivalent to 1)44(1 provided we identify QtZ — AttKax. The analogy between 
the extremal points of </> and the Coulomb-gas is not perfect and is valid only in the 
high-temperature region (see |29l 15 Ip. As already pointed out in section ITTl the 
scalar field is a nonsingular, single-valued object and must be distinguished from a 
field of polar angles 9. The latter field (not considered here) is multivalued and has 
singularities which behave exactly as Coulomb charges j521 [531 . 

For a flat geometry Gij = Sij, the average defect density is zero, since 

is also a flat metric with vanishing curvature. 



4. The two-point function 

The calculation of the two-point function C(r) = (p(r)p(O)) is done here for a 
Euclidean d-dimensional geometry and a particular class of (isotropic) vector fields 
Ui with independent components ui and U2. The correlation function reads 

Mr)u,{r'))=S,,G{\r-r'\) (48) 

equipped with an appropriate short-distance cutoff. Examples are random waves [SJ 
with G{r) = Jo{kr) and the Gaussian XY-model with G{r) = KQ{y/Tr), where 
Jo is zeroth order Bessel function and Kq is the modified Bessel function of order 
zero. Gradients of Gaussian scalar fields do not fall into the above class, since the 
components of the gradient di4> are not independent. The metric tensor and its inverse 
are now 

9iaJI3 = Sij {SapGo + CapG{r)) 
rc.,jl3 _ e. . ^cpGo - Co,pG{r) 

where Gq = ui + G'(O), Cap = 1 — 5ap, i,j = 1 . • . d and a,(3 — 1, 2. The determinant 
of g is det g — (Gq — G^(r))'^. The affine connection is 

^k'y,ia;mtL = ^07 {diUk{r^)um{r ^1,)) = 5a^e^^i5kmdiG{r) (50) 

where 7 is not summed and ei2 = 1, 621 = —1, en = 622 = 0. We obtain the curvature 
tensor 

Rki,l\,ia,jP = Ski {SajSpx — Sa\Spj) {~~didjG{r) 

-e^^exud,G{r)d,G{r) ^2 _ g2(^) 
= - SkiCape^x I didjG{r) + diG{r)djG{r) 



G{r) 



Gl-G^{r) 

^Gl- G2(r) 6Mec.pe^xd,d,K{r) (51) 
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where K{r) = arcsin(G(r)/Go) and the indices 7, A are not summed. The two-point 
correlation function C(r) = {p{r)p{Q)) now reads 

{2^YC{r) = {Gl - G2(r))-^/2 / {di^AOf'^ 



X exp [el6U7f^^Gl-G^{r)d.d,K{r) 

= {-lYd\ J (dr;2dr;i)^exp {rjlrjp,d,Kir)) 

^{-lfdldet{^^^jK{r)) (52) 

since the term (. . .)'^/d\ of the exponential is the only one which contributes to the 
result. The square root -^/Gq — G^(r) is therefore cancelled by the prefactor. The 
integration over the 6'-variables yields an extra factor {—l)'^dl. Denoting d{. . .)/dr = 
(. . .)' we end up with 

{2nYG{r) = {-lYd\K{r)" {K{r)' /rf'^ . (53) 

This result was already obtained by direct calculations ^| I28L I4UI [5^ for the limiting 
case w ^ 0. The correlation function G(r) is a functional of the ratio G(r)/Go and 
is therefore independent of the scale T of the field Ui (in the limit of a sharp delta 
function w — > 0). 

^.1. Singular behaviour of correlations for coinciding points and neutrality 

Each vector field placed on a closed surface is subject to the topological constraint 
l[T]l. implying a neutrality condition for correlation functions 

/ dV^/d^ (p(fT)p(a')) - {p{<y') J2 9") = 2(1 - 7) (P(^')> , (54) 

where Gij (a) is the metric tensor of the surface. Analogue constraints apply to higher 
order correlation functions. To find a neutrality condition for vector fields on an 
infinite planar geometry, we note that the density p can be written as a divergence 
(for a finite width w) 



-■ii,...,id 

k 



Yl d,,Uk {2ttu;)-^/^ exp {-{ukf / {2lj)) 

k 

J|5,,erf(Mfe/Vw) 



— <-ii,...,id 

fc 

= di ^ei,i2,...,i^erf('Ui/Vt^) di^eTi{uk/ V^)^ (55) 

where erf(a;) = (27r)~^/^ Jp" ds exp(— s^/2) is the error function and eij_....i^ is the 
antisymmetric unit tensor in d dimensions. Then the total charge of a vector field is 
zero, i.e. 

dV {p{r)p{r')) = (56) 

provided that the correlations (erf(Mi(0))erf(uj(r))) ~ {ui{0)uj{r)) decay fast enough 
(at least ~ r"^/^) for r 00. On the other hand, the correlation function should 
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behave like a delta-function for r ^ r' and — > since 

{P{r)p{r')) = J2 - IP^'i^' - ^0)) 

= 6'{r - r') ^ {5\r - r„)) + ^ (g,<52(r - r„) qp5\r' - r^)) 

where no = (| Aei{diUj)\5'^ [u)) is the absolute density of zeros We consider now a 
planar geometry and vector fields with an isotropic distribution. Then rig is spatially 
constant and both {p{r)p{r')) and '}Za±B {Qa^^if ~ ^q) 9/3^^(^' ~ '"z?)) are functions 
of |r — r'l only. We integrate equation (|57|) over r' and obtain (see also j28l 1551 1551 14(jj ') 

= no+ / dVC(r) (58) 

J \r\>a 

where we have used equation (|56|l . The integral extends over two-dimensional space 
with the exception of an infinitesimal small disk (radius a — > 0) at the origin. For 
d = 2 and vector fields with correlation function (|48|l we have found (see equation (|53|l 
in the preceding section) 

{(K'(r)f\' 
(27r)^r 

The correlation function of the vector field u must be regularized at short distances and 
can therefore be expanded for small radu G(r) = G(0) -I- r^G" (0)/2 . . . (G"(0) < 0), 
which implies (for small w) 

K{r) = arcsin(G(r)/Go) = f ^ (|j + ' 



TT /-G"(0) 



(60) 



2 V ^0 

where = 2w/(-G"(0)). We obtain for w -> with K'{r) ^ for r ^ 00 



1 



-no = 2TT rdrC{r) = — {K' {r)Y - TTT^.- (61) 



27r 



00 



G"(0) 



0+ 27rG(0)' 



The function K{r) has for — > a conical singularity at the origin. The correlation 
function G(r) cx A.et{didjK{r)) can now be interpreted as the Gaussian curvature of 
that cone, which is in fact concentrated at its tip. Using equation 159(1 . we find the 
anticipated short-distance behaviour of the two-point correlation function G(r) 

^, , (~G"{Q)\ 2a2 -G"(0).2, , .2, ^ 
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5. Conclusions 

We have proposed a scheme to evaluate correlation functions of the (charge-) density 
of zeros p of random vector fields Ui with a Gaussian distribution. The zeros carry 
a charge q = sign J = ±1, where J = det (diUj) is the Jacobian of the vector field 
right at the zero. In two dimensions for instance, positive zeros resemble sources, 
sinks or vortices, while negative ones have a saddle-like flow. We could show that the 
correlation functions {p{ri) . . .p{rf)) can be expressed through the curvature of an 
abstract / x ci-dimensional Riemann-Cartan manifold (d is the dimension of both the 
vector field and the embedding space). In the case of gradient fields Ui = di^, the 
correlation functions are the total curvature of a certain / x rf-dimcnsional Riemannian 
manifold. As an application, we have calculated the mean charge density for various 
vector fields on curved two-dimensional surfaces. Furthermore, we derived the two- 
point function for isotropic vector fields with independent components. We checked the 
neutrality sum rule, which follows from the topological constraint of the total charge 
and obtained in addition the short-distance behaviour of the two-point correlation 
function. We have seen that the zeros of Gaussian vector fields are a model for the 
two-dimensional Coulomb gas in the high-temperature region. This opens a convenient 
route to study these systems in curved or bounded geometries. The geometrical view 
helps to calculate one- and two-point correlation functions and could also be rather 
useful to 'tame' the higher order correlation functions of Gaussian zeros. 

The present work demonstrates that the signed zeros of vector fields with a 
Gaussian distribution are not only sensitive to curvature, the correlation functions of 
the charge density are the curvature of an abstract manifold, which is closely related 
to the embedding manifold. 
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